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ABSTRACT 

We report the discovery of a young and energetic pulsar in the Parkes multibeam survey of the Galactic plane. 
PSR J1016-5857 has a rotation period of 107 ms and period derivative of 8.0 x 10"'"*, implying a characteristic 
age of 21kyr and spin-down luminosity of 2.6 x lO-'^ergs"'. The pulsar is located just outside, and possibly 
interacting with, the shell supernova remnant G284.3-1.8. Archival X-ray data show a source near the pulsar 
position which is consistent with emission from a pulsar wind nebula. The pulsar is also located inside the error 
box of the unidentified EGRET source 3EG J1013-5915, for which it represents a plausible counterpart. 

Subject headings: ISM: individual (G284.3-1.8, 3EG J1013-5915) — pulsars: individual (PSR J1016-5857) — 
supernova remnants 



1. INTRODUCTION 

Observations of young rotation-powered pulsars offer unique 
opportunities to illuminate several fundamental questions about 
the physics of neutron stars and supernova remnants (SNRs). 
In particular, synchrotron emission from a pulsar wind nebula 
(PWN) can be used to determine properties of the relativistic 
pulsar wind, to probe the density of the ambient medium, and to 
study the evolution of SNR ejecta and their interaction with the 
local ISM. Confirmed associations between pulsars and SNRs 
also provide the opportunity to directly constrain the birth prop- 
erties of neutron stars, including their space velocities and ini- 
tial rotation periods. Such associations are often difficult to 
judge on the basis of the distance and age estimates for both 
objects, and a very powerful method of confirming an associa- 
tion is the observation of interaction between pulsar/PWN and 
SNR. However, young pulsars are rare, and cases where such 
interaction is apparent are rarer still. 

One of the m ain aims of the Parkes rn ultibeam survey of the 
Galactic plane ( Manchester et al. 2001 ) is precisely to detect 
young pulsars, and at least ~ 30 of the 600 pulsars discovered 
to date are youthful, with characteristic age r,. < lOOkyr S ev- 



eral of these are energetic as well (e.g., D'Amico et al. 2001 ). 
In this Letter we report on one such object, which archival ra- 
dio. X-ray, and 7-ray observations suggest is remarkable. 

2. RADIO PULSAR J10I6-5857 

PSR J1016-5857 was discovered in data acquired on 1999 
Marc h 15, collected in the cour se of the pulsar multibeam sur- 
vey ( Manchester et al. 2001 ) at the 64-m Parkes telescope 
in NSW, Australia. The pulsar has been the subject of regu- 
lar timing observations at Parkes since 1999 May. These use 
the center beam of the multibeam receiver and the same data- 



acquisition system as the survey to collect total-power samples 
every 0.25 ms from 96 contiguous frequency channels between 
1230 and 1518 MHz. The 1-bit digitized samples are written to 
tape for subsequent analysis, together with relevant telescope 
parameters, including the start time of the observation referred 
to the observatory time standard. The time on-source is 9 min 
on average each observing day. 

The timing data are analyzed in standard fashion: samples 
from higher frequency channels are appropriately delayed to 
compensate for dispersion in the ionized interstellar medium. 
Data from all frequency channels are then summed and the re- 
sulting time series is folded modulo the predicted topocentric 
pulsar period to generate an average profile for each observa- 
tion. A sum of several aligned average profiles yields the "stan- 
dard profile" (Fig. |l]), with which individual profiles are corre- 
lated to obtain times-of-arrival (TOAs). 

The TOAs are analyzed with the TEMPO software package^: 
they are converted to the solar- system barycen ter using the 
JPL DE200 planetary ephemeris ( |Standish 199C| ), and fitted to 
a Taylor series expansion of pulse phase by minimizing the 
weighted rms of timing residuals (observed minus computed 
phase). This yields precise values of pulsar rotation frequency 
u, its derivative i), and celestial coordinates. Many young pul- 
sars, including PSR J1016-5857, experience rotational insta- 
bilities, and in the presence of such "timing noise" particular 
care must be taken in order to minimize contamination of fit- 
ted parameters. Following the prescription of Arzoumanian et 
al. (1994) we first "whiten" the timing residuals, by fitting for 
position, i^, i>, and higher frequency derivatives as required (z>' 
in this case). The resulting celestial coordinates and uncertain- 
ties are our best estimates of these parameters and are listed in 
Table |l| as are v and 1), which we determine from one addi- 
tional fit with position held fixed and j>' set to zero. We then 
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perform one final fit which includes iJ. This is not a term de- 
scribing the secular evolution of the neutron star rotation, but 
results from rotational instability in the pulsar. In fact, both the 
sign and magnitude of i) strongly suggest that it mostly repre- 
sents recovery fr om a glitch that occurred prior to the discovery 
of the pulsar (cf. |Lyne 1996 ). The timing residuals correspond- 
ing to the solution presented in Table |l], based on 1 .5 yr of data, 
are shown in Figure [l|. 
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Fig. 1. — Pulse profile and post-fit timing residuals for PSR J1016-5857. 
CoiTesponding parameters and weighted rms JEsiduals for "red" (top) and 
"white" (bottom) solutions are listed in Table hi Top: (inset) average pulse 
profile at 1400 MHz, with 0.6 ms resolution. The full profile width at half the 
maximum intensity is »'5o = 3.3 ms, and at 10% it is n'lo = 11 ms; (main) fit for 
u and v only (position held fixed at best fitted value). The quasi-cubic trend in 
residuals indicates presence of timing noise. Bottom: as for upper panel, with 
additional fit for v (see § 0). 

With period P = \/v = 107ms and a high P, PSR 
J1016-5857's inferred characteristic age, r^ = P/2P = 21 kyr, is 
relatively small and its spin-down luminosity, E = Att^IP/P^ = 
2.6 X lO^'^ergs"', large. The surface dipole magnetic field 
strength is B = 3.2 x W^iPpyl^ = 3.0 x 10'^ G (Table |I]). These 
parameters make it similar to the small number of "Vela-like" 
pulsars, which have r^ ^ 10"^-10^ yr and E > 10^^ erg s"' . Prior 
to the Parkes multibeam survey, only nine such pulsars were 
known. Thus far we have already doubled this number'^. Pul- 
sars in this group have high levels of glitch activity, and many 
are associated with PWNe and/or SNRs. 

3. SUPERNOVA REMNANT G284.3-1.8 

PSR J1016-5857 is located, at least in projection, near the 
supe rnova remnant G284.3-1.8 (MSH 10-53; M ilne et al. 
1989; see Fig.g). SNR G284.3-1.8 is an incomplete radio shell 
with non-thermal spectrum and wi th significant polari zation, 
interacting with molecular clouds (Ruiz & May 1986). The 



pulsar is located exactly at the very tip of a "finger" of emission 
on the western edge of the remnant, and ^ 15' from the SNR's 
approximate geometrical center. 

In the Galactic plane, the spatial density of SNRs and pul- 
sars is high, so the fact that a young pulsar lies near the edge 
of a SNR does not necessarily imply that they are associated: 

'See http://www.atnf.csiro.au/research/pulsar/pmsur\ . 



the probability of finding a pulsar lying by chance within 15' 
of a SNR center i n this region of the Galaxy (270° < I < 300°) 
is about 5% (cf. Green 2000). However, the morphology of 



this particular pairing appears to bear a resemblance to that 
of PSR B 1757-24 (another Vela-like pulsar) and its associated 
SNR G5.4-1.2 ("the Duck"). In the latter case, the pulsar has 
apparently caught up with, penetrated and overtaken the rim of 
its SNR, leaving a trail in its wake; the pulsar is now embed- 
ded in a ra m-pressure confined (bow-shock) PWN (M anchester 
et al. 1991; [Frail & Kulkarni 1991[ ). Because of the Umited 
resolution and sensitivity of the MOST data shown in Figure g, 
it cannot yet be determined whether a similar association exists 
between PSR J1016-5857 and SNR G284.3-1.8. Nevertheless, 
the location of PSR J1016-5857 at the tip of a bright structure 
connected to G284.3-1.8 is suggestive of a physical association 
between the pulsar and SNR. 
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Fig. I? — MOST radio map at S43MH7 nf SNR mS4 3-1 S (.M ilne et 
al. 1989), with a spatial resolution ot 43" . jJie position ot FSK J 1016-5857, 
known with sub-arcsecond accuracy (Table Mj), is indicated with a large cross 
at the tip of the protuberance on the western edge of the SNR. 



If the pulsar and SNR are associated, the pulsar is at the dis- 
tance of the SNR, c/ = 3 kpc, on the near side of the Carina spiral 
arm, as determined from optical and CO line observ ations with 
an uncertainty of only « 20% ( Ruiz & May 198q ). This dis- 
tance is smaller than that inferred from the pulsar's dispersion 
measure (DM) and a mo del for the free electron distribution in 
the Galaxy, d = 9;!;2kpc (Taylor & Cordes 1993). However, the 
model is at its most uncertain along spiral arm tangents, the di- 
rection of PSR J1016-5857, so that the discrepancy may not be 
as serious as it appears. A similar example is provided by the re- 
cently d iscovered ^ 1600-yr-old pulsar Jl 1 19-6127 (C amilo et 
al. 2000) a nd its clearly associated SNR G292.2-0.5 ( Crawford 
et al. 2001; pivovaroff et al. 2001 ), where the model puts the 
pulsar atd> 30 kpc, when in fact it lies at t/ <; 8 kpc. Hypoth- 
esizing the PSR J1016-5857/SNRG284.3-1.8 association, we 
consider d = 3 kpc hereafter 

The age of G284.3-1.8 is estimated by Ruiz & May (1986) 
as r^ 10 kyr, by fitting line profiles observed from an opti- 
cal filament to models of interstellar radiative shocks. This 
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is similar to the pulsar's characteristic age, r^ = 21 kyr. How- 
ever, pulsar characteristic ages are not necessarily accurate es- 
timates of true age t: e.g., in the B17 57-24/G5.4-1.2 sys- 
tem, T ^ Tp ( jaensler & Frail 2000[), while for t he pul- 
sar in SNR Gll.2-0.3, t <s: t^ (|Kaspi et al. 200ll). Still, 



if the PSR J1016-5857/G284.3-1.8 association is real, and 
the age of the system is ^ r^, we can estimate the veloc- 
ity required for the pulsar to move ballistically from its pu- 
tative birthplace (SNR center) to its present location: v ^ 
500((i/3kpc)(21kyr/T)kms"'. This velocity is within th e 
range expected for young neutron stars ( Ly ne & Lorimer 1 994 ) , 
but is highly supersonic in the local ISM. A characteristic bow- 
shock PWN is then expec ted, as is observed for PSR B 1757-24 
(Gaensler & Frail 2000). Such a bow-shock nebula, of typ- 
ical extent <; 1', could not be seen in Figure g. Further- 
more, the larger-scale finger connecting to the SNR need not 
point in the direction of motion, since the "wake" left be- 
hind by the pulsar can be distorted by ISM motions and turbu- 
lence. Finally, in comparing the PSR J1016-5857/G284.3-1.8 
and PSR B1757-24/G5.4-1.2 systems, we note that the pul- 
sars are uncannily alike; P = 107/128 ms, r^ = 21/16kyr, E = 
2.6/2.5 X 10^''ergs-',£ = 3/4x lO'^G. 

4. X-RAY OBSERVATIONS 

Energetic young pulsars such as PSR J1016-5857 often pro- 
duce detectable high-energy emission, especially if their dis- 
tance is not too large (as we have argued in the previous sec- 
tion). A search of the HEASARC archive for serendipitous 
high-energy observations of PSR J1016-5857 has yielded one 
pointing at SNR G284.3-1.8, obtained with the Imaging Pro- 
portional Counter (IPC) on-board the Einstein Observatory. 
The IPC is sensitive to 0.1^.5keV X-rays over its 1° x 1° 
field-of-view with an image resolution of l'-3'. The standard 
processed image from this 1.9ks exposure (Seq. 15288), ac- 
quired on 1980 February 5, is shown in Figure H. The pul- 
sar position is near the center of the IPC field and coincides 
with an excess of localized X-ray emission, the most signifi- 
cant in the field. This emission is consistent, given the limited 
number of counts, with a point source, with possible extended 
emission. From the original event data files, we extracted 134 
counts from a r = 3.'0 circle around the peak of emission and 
estimated the background counts within a surrounding annulus 
at 3.'2 < r < 6f2. The 45 background-subtracted counts yield a 
source significance of 4.8 a and a centroid which is -^ 1 .'8 offset 
from the radio position. While much larger than the nominal 
Einstein pointing error (<; 10"), this offset is consistent with 
the IPC instru mental measureme nt error and the small number 
statistics (e.g., Harris et al. 1990| ). 

To estimate the luminosity of the source we assume a typi- 
cal pulsar power-law model with photon index w 2, as might 
be appropriate for magnetospheric emission or origin in a 
non-thermal nebula. A measure of the Galactic neutral hy- 
drogen column density to the pulsar is estimated from the 
Hi study of Dickey & Lockman (1990)'". This gives Nh « 
10^^ cm"^, consistent with the rule of thumb of 0.1 free elec- 
trons per H atom and the pulsar's DM. For isotropic emission, 
we obtain an unabsorbed X-ray luminosity of Lr,o. 1-4.5 kev ^ 
10^^(t//3kpc)^ergs"'. This represents an efficiency for con- 



version of rotation-driven luminosity into X-rays of L^/E ~ 
5 X 10"^ in the Einstein band. Such an efficiency for the point 
source and any compact nebula surrounding it is consistent 
within uncertainties w ith that observed for young pulsars (e.g., 
Seward & Wang 1988|). 
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Fig. 3.— Einstein IPC observation of the PSR J1016-5857/SNR 
G284.3— 1.8 field: a possibly exteiided source is detected near the position 
of the pulsar (cross-hairs; see § W). The ellipse overlaid represents the 
approximate la mnmnr nf positional uncertainty for the EGRET source 
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Given the limitations of the detection, the origin of the ob- 
served emission is not certain at present. Possible contributions 
include emission from: a point source (e.g., pulsed magneto- 
spheric emission); a compact nebula (bow-shock or static); the 
SNR; or even interaction of the SNR with a molecular cloud (cf . 



Ruiz & May 1986)". However, the inferred flux is in excess of 
the thermal emission expected from a '--^ 20 kyr-old neutron star 
for a reasonable Nh, while the measured flux is approximately 
as expected from a PWN in such a pulsar system. Also, if the 
extended emission is real, it appears to point back towards the 
center of the SNR, suggestive of a ram-pressure confined wind 
nebula. 

5. 7-RAY SOURCE 3EGJ1013-5915 

As shown in Figure |, PSR J1016-5857 and the Einstein X- 
ray source are within the « 43' error box of EGRET source 
3EG J1013-5915. This 7-ray source has flux (3.3 ±0.6) x 10"^ 
photo ns cm~^s~' (> 100 Me V ), and photon index T = 2.32 ± 



0.13 (Hartman et al. 1999). This photon index is slightly 



larger than for most 7-ray pulsars, but consistent with those 
of the Crab (T = 2.19), and the Vela-like PSR J2229H-6114 
(F = 2.24), which is al most certainly the so urce of the 7-ray s 
from 3EG J2227H-6122 ( |Halpernet airiOOlb. 3EG J1 013-5915 
is possibly a composite source ( [Hartman et al. 1999 ), in which 
case the photon index of any pulsar component may be differ- 



Accessible via the nH tool at ittp://heasarc. gsfc.nasa.gov/docs/corp/tools. html 



' ' Only one cataloged object exists within 2' of the pulsar position, the A2 star HD 
while if it is a main-sequence A2V star it lies at < 1 kpc. At such a distance, the flux 
known for any A2V star. 



302636, 1' away. If this star were a supergiant it would lie outside the Galaxy, 
measured from the Einstein source would imply a luminosity in excess of that 
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ent from the nominal value, and its flux will be lower than that 
liste d. The source is not obviously variable (M cLaughlin et al. 
1996), and the flux integrated from lOOMeV to lOGeV with 
the nominal photon index corresponds to an isotropic luminos- 
ity L^ - 1.6 X 10^^V/3kpc)2ergs-i,or-0.06£'. 

The chance probability for superposition of a Vela-like pul- 
sar and an EGRET error box in the Galactic plane is not negli- 
gible: in the 1500 square degrees covered by the Parkes multi- 
beam survey, where there are about 10 Vela-like pulsars known, 
EGRET error b oxes cover 24 square degrees or 1.6% ( Hart- 
man et al. 1999). The chance probability that at least one of 
those pulsars is aligned with one of the error boxes is there- 
fore 1-(1-0.016)'" = 0.15. The actual probability is somewhat 
smaller if one considers only non-variable EGRET sources, but 
it is nevertheless significant. 

However, both L^/E ^ 0.06 and L^/Lx ~ 100 are ap- 
proximately as expected for a pulsar of the age and E of 
PSR J1016-5857 ( e.g., [Thompson et al. T999|, and Se- 
ward & Wang 1988), as is the lack of variability in 7-rays. 
Also, every pulsar with a value of E/cfi greater than that of 
PSR J1016-5857 has been detected in 7-rays (although the dis- 
tances of PSR J1016-5857 as well as those of other pulsars are 
substantially uncertain). Altogether, this makes an association 
between PSR J1016-5857 and 3EG J1013-5915 possible and 
intriguing. The proof that 3EG J1013-5915 is the 7-ray coun- 
terpart of PSR J1016-5857 would come from the detection of 
7-ray pulsations. This will be difficult with EGRET data, since 
the timing noise of the pulsar and the likelihood of glitches (§ 0) 
prevent the reliable backwards extrapolation of the ephemeris 
required to search for 7-ray pulsations. 

6. DISCUSSION AND FUTURE WORK 

We have discovered an energetic and youthful pulsar, ap- 
parently located just outside SNR G284.3-1.8. Based on 



its spin characteristics, PSR J1016-5857 belongs to the rare 
class of Vela-like pulsars, those with r^ ~ 10^-10^ yr and E > 
10^^ ergs"'. These pulsars are of particular interest, as many 
are embedded in PWNe and are observed to be interacting with 
their SNRs. Additionally, such pulsars are very efficient at con- 
verting their spin-down luminosity into 7-rays, and form the 
majority of EGRET sources identified with pulsars. 

We have argued, based on its positional coincidence with 
one notable feature in the SNR, that PSR J1016-5857 may 
be interacting with G284.3-1.8, and therefore that the two ob- 
jects are physically associated. Further supporting this notion 
is the detection of an X-ray source nearby PSR J1016-5857 
and the positionally coincident unidentified 7-ray source 
3EG J 101 3-59 15. The luminosities of both X- and 7-ray 
sources, if located at the distance of G284.3-1.8, are consis- 
tent with provenance from PSR J1016-5857. 

Confirming this picture will require high-resolution radio and 
X-ray observations that show the proposed interaction between 
pulsar and SNR. If a bow-shock is detected, its detailed study 
may constrain the pulsar velocity, age, spin-evolution, and am - 
bient medium (cf. [Gaensler & Frail 20()0t [Olbert et al. 200l|). 



Confirmation of the nature of the EGRET source may have to 
await future missions such as GLAST. 
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Table 1 
Parameters of PSR J1016-5857 



Parameter 



Value 



R.A. (J2000) 10''16'"2in6(l) 

Decl. (J2000) -58°57'12;'1(1) 

Rotation frequency, v (s"') 9.31274229245(4) 

Period, P (s) 0.1073797565311(5) 

Frequency derivative, (s"-) -6.991787(7) X 10"'^ 

Period derivative, P 8.061818(8) X 10"'" 

Second frequency derivative, i) (s"^)1 1.03(2) X 10"^^ 

Epoch (MJD) 51730.0 

Dispersion measure, DM (cm"' pc) . . 394.2(2) 

Data span (MJD) 51451-52002 

Rms timing residual (red/white) (ms)*" 4.3/0.24 

Flux density at 1400 MHz (mjy) .... 0.46(5) 
Derived parameters: 

Galactic longitude, / (deg) 284.08 

Galactic latitude, b (deg) -1.88 

Distance, d (kpc)*-' 3? 

Characteristic age, Tc (yr) 2.1 X 10* 

Spin-down luminosity, E (erg s"') . 2.6 X 10^* 

Magnetic field strength, B (G) 3.0 X lO'^ 



Note. — Numbers in parentheses represent 1 u uncertainties in 
the least-significant digits quoted. 

"fitted for separately (not a secular term) — see § Q. 

''see Figure ry. 



